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Abstract. We consider the statistics of volume fluctuations in a one-dimensional 
. classical gas of non-interacting particles confined by a piston, and subjected to an 

arbitrary external potential. We show that despite the absence of interactions between 
particles, volume fluctuations of the gas are non-Gaussian, and are described by 
generalized extreme value distributions. The continuous shape parameter of these 
distributions is related to the ratio between the force acting on the piston, and the force 
g ■ acting on the particles. Gaussian fluctuations are recovered in the strong compression 

O \ limit, when the effect of the external potential becomes negligible. Consequences for 

the thermodynamics are also discussed. 

<N : 
> ' 

as : 

y—( ', 1. Introduction 

; 

■ Extreme value distributions, describing the fluctuations of the k th largest value in a 

■ set of random variables [U [2], have recently been shown to play a role in physics 
well-beyond their standard area of application. Such distributions are indeed found 
analytically in different problems where extreme values are not involved a priori, such 

^ ■ as global measures in 1 / /-noise P, H] , in interfaces in random environment [5] , in coupled 
o5 ' quantum oscillators [6], or in the level density of Bose gases [7j and related problems 
of integer partitions [8]. Moreover, distribution functions describing the fluctuations of 
disparate global measures in complex systems [9] can be well approximated by extreme 
value distributions generalized to real values of the parameter k [TUl [TT1 [T2] . Examples 
also include Burgers turbulence [13], models of heterogeneous glassy dynamics [HI [15] . 
relaxing granular gases [16], and very recently, order parameter fluctuations in a liquid 
crystal close to a critical point p2]. Generalized extreme value distributions were also 
found analytically in a stochastic cascade model with dissipation |18j . 

A scenario has been proposed recently to explain the rather puzzling emergence 
of extreme value distributions in contexts not clearly related to extreme processes 
[T8l [19] . One can actually reformulate the original problem of extremes as a problem 
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of sums of non-identically and generically correlated random variables, leading to the 
same extreme value distributions. Provided that the joint probability of the summed 
random variables has a particular form, extreme values could result from random sum 
problems. For instance the appearance of a Gumbel distribution in the l//-noise model 
is understood in a simple way within this framework [18]. Furthermore this mapping 
provides a new view on generalized extreme value distributions. It is indeed possible 
to generalize the equivalent sum problem simply by extending the joint distribution to 
real values of the index k: generalized extreme value distributions then find a natural 
interpretation as limit distributions of sums of non-identically and generically correlated 
random variables [HI [19], with a particular form of joint probability (see Eq. (TlOl) below). 
However, with the exception of uncorrelated but non-identically distributed variables 
which further generalizes the 1/ /-noise model by including a low frequency cut-off [19|, 
the aforementioned class of correlated random variables looks rather formal and far from 
physical applications. 

In this paper, we illustrate on a very simple statistical model, namely a one- 
dimensional classical gas of independent particles, how the above class of random 
variables can find a natural application in a physical context. This type of model 
has a long history and has proved very useful for the illustration of statistical concepts 
[20| |2T| [22l [23l [2U [25]. However, at variance with most previous studies on such a 
model, we consider here a gas confined both by a piston and by an external potential. 
We find that although the system is composed of independent and identical particles, 
the distributions of volume fluctuations are non-Gaussian, and described in the limit of 
a large number of particles by generalized extreme value distributions, indicating that 
correlations appear in the syste The nature of the correlations is however rather 
subtle: the local volume between particles is correlated, but their positions are not. 

Note that connections between the one- dimensional classical gas of particles (or 
Jepsen gas) and extreme value statistics have already been reported when the gas can 
freely expand, in the absence of piston and potential [25] . However, these results concern 
the velocity of the rightmost particle rather than the volume of the gas, and are thus of 
a different nature. 

2. A simple model of a one-dimensional confined gas 

We consider an ideal system composed of N point-particles placed in a cylindrical 
container with long axis z and with a small diameter with respect to its length (quasi 
one-dimensional geometry). The model is illustrated in a schematic way in Fig. [IJ The 
position of the particles along the z-axis is denoted as Zi, % — 1, . . . , N. A hard wall, 
acting as a reflecting boundary is placed at z = 0, so that particles are constrained to 
remain on the half-space z > 0. The container is closed by a piston that can freely 
move along the z-axis. The position of the bottom of the piston is denoted z p , so that 

| Another reason for non-Gaussian fluctuations would be that the variance diverges, leading to Levy- 
stable laws [25], but this is not the case here. 
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Figure 1. Schematic illustration of the model: point-particles at position z% are 
placed in a container closed by a moving piston, at position z p . An external potential 
(illustrated as the gravity g ) acts on the particles, and a force f p (that may be of the 
same origin as that acting on the particles) is applied on the piston. 



the volume of the system is given by V = Sz p , where S is the cross-section of the 
cylinder. In addition, particles are subjected to an external potential U(z), z > 0, 
that tends to confine them in the small z region, that is U(z) is assumed to be an 
increasing function of z. The piston is also subjected to an external potential U p (z p ), 
that may differ from the potential acting on the particle. The potential U p may for 
instance be the gravitational potential acting on the piston. It may also be the potential 
caused by an operator exerting a constant force f p on the piston, corresponding to a 
linear potential U p (z p ) = —f p z p , or by a spring fixed on the piston, corresponding to a 
quadratic potential U p (z p ) = \K(z p — z ) 2 - U p (z p ) may also result from the superposition 
of the different types of potentials mentioned above. 

The container plays the role of a heat bath that thermalizes the particles at a given 
temperature T. Collisions between the particles and the piston are elastic, so that the 
piston in turn thermalizes at temperature T (the piston does not have any internal 
structure, only its translation degree of freedom thermalizes). Note that the present 
model differs from the so-called Jepsen gas [201 ED E21 E31 El] , due to the presence of a 
heat reservoir and of an external potential. 

3. Distribution of volume fluctuations 

We now characterize quantitatively the volume fluctuations of the system through the 
position, z p , of the piston. When all the particles and the piston are equilibrated at 
temperature T, the equilibrium distribution reads, with z iy z p > 0, 

1 N 

P N {z u ...,z N ,z p ) = - e-^pfeO J] e~ mZi) ©(^p ~ *), (!) 
L i=i 

where (3 = 1/ksT is the inverse temperature, Z is the partition function and O(x) is 
the Heaviside function. The volume distribution Qn( z p ) is obtained by integrating over 
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the variables Z\, . . . , z N , leading to 

Qn(z p ) = le-^-W (frdze-Wy. (2) 

Before going into more detailed calculations, we briefly discuss two simple limiting cases 
of interest, namely U(z) = and U(z) = U p (z). The case U(z) = is the most standard 
equilibrium case, for which 

Q N (z p ) = ^e-^\ (3) 



and Gaussian fluctuations are recovered in the large N limit 

In contrast, if the potential acting on the piston is the same as that acting on the 
particles, one can express P(z p ) in the following form: 

Qn(z p ) = ^G(z p ) N+ \ (4) 



dz p 



with 



Jo p dz e- pu{ A 

G{Zp) = s™dze-m*r (5) 

Volume fluctuations then exactly map onto an auxiliary problem of extreme values, 
namely, the fluctuations of the maximal height of N + 1 independent particles (with 
no piston) with positions z[ in a potential U(z). The function G(z p ) is simply the 
probability that the position of a particle subjected to the potential U(z) is smaller 
than z p . Hence, G(z p ) N+1 is the probability that the positions of N + 1 particles are less 
than Zp, which is nothing but the cumulative distribution of the maximum of the N+l 
positions z\. From Eq. (TJJ, it follows that P(z p ) is the distribution of max(z[, . . . , z' N+1 ), 
so that in the large N limit, fluctuations of z p are described by standard extreme value 
distributions. 

In this paper, we are mostly concerned with the intermediate situation where 
particles are submitted to a force derived from a potential, but where this force is 
smaller than that acting on the piston. A typical situation of this type is that of a 
system placed in a gravitational field, as the piston generically has a mass larger than 
that of the particles. 

Let us now compute the asymptotic volume distribution for general potentials U (z) 
and U p (z). One way to tackle this issue could be to start directly from Eq. (j2J). Rather, 
as a short-cut, we take an alternative approach using the results derived in [19]. To 
this aim, we introduce the intervals between the ordered positions of the particles in the 
following way. For a given set of values z\, . . . ,z^ satisfying < Zi < z p for all i, we 
introduce a permutation a of the integers 1, . . . , N such that z a (i) < z a (2) < • • • < Za(N), 
and we define the space interval hi between particles through 

hi = z a{i) - z ff (i-i), i = 2, . . . , N. (6) 

For convenience, we also introduce the variables h± and h^+i, 

hi = Zo-(i), h N+ i = z p — z a (jsi). (7) 
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It is then straightforward to express Zi as a function of the variables hj, namely 

o- -1 (*) N+l 

z i= E h v (* = 1> ■ ■ ■ > z p= h r ( 8 ) 

i=i 3=1 

where a~ l is the inverse permutation of a. The system can then be described by 
the set (h±, . . . , /ijv+i), with hi > for all i, up to an arbitrary permutation of 
the N distinguishable particles. A given set (hi, . . . , /ijv+i) then corresponds to AM 
configurations of the particles (the position of the piston is fixed when the h^s are 
given). Summing over the corresponding AM configurations (zi, . . . , z^) in Eq. (CQ), one 
obtains the equilibrium probability distribution 

N 



P N (hi 



h 



N+l/ 



M e -Pu P (E 



X > h - \[< 
k=l 

-li 



(9) 



where we relabelled the factors in the product using k = a (i). Note that this last 
equation is also obtained in the case when particles cannot cross each other. In this case, 
the 'no-crossing' constraint needs to be taken into account from the outset in Eq. (Tj[|), 
and Eq. (Q is rather straightforwardly obtained, since there is no need for reordering the 
positions of the particles. Note also that the case of indistinguishable particles leads to 
a result similar to that of distinguishable particles that cannot cross each other, since in 
both cases it is not possible to generate a different configuration through a permutation 
of the particles. 

The distribution P^(hi, . . . , h^ + i) given in Eq. turns out to be quite similar to 
the joint distribution describing the class of correlated random variables introduced in 
[T9] . Up to slight notation changes [jj, the latter reads: 



Jtv(/ii, /iv+i) — Km £1 




N+l 



(10) 



1=1 



u- 1 



Q(F) is an arbitrary function of F, with < F < 1, and the function F(z) is defined as 



F(z) 



P(z')dz'. 



(11) 



The function P(z) has the properties of a one- variable probability distribution, namely 
it is a positive function such that P(z')dz' = 1. In order to map the gas model onto 
Eq. (flUj) . we make the following identification: 

P(z) = Xe- (3Uiz \ (12) 



n[F(z p )}P(z, 



(13) 



for all z, z p > 0, and where A is a normalization factor. With this identification, Eq. 
can be rewritten as 



P N (h 



i- 



h 



N+l, 



AM 



n 




N+l 



Up E*i 



(14) 



k=l 



U'=l 



§ Starting from the distribution J/v(ui, . . . , un) defined in |19j . we change N into N+l, and reverse 
the order of index, that is, we define the variables hi = ujv+2-i- 
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which is precisely the same form as in Eq. (1101) . As P(z) > for all z > 0, it results 
from Eq. f[TT|) that F(z) is a strictly decreasing function of z > 0, so that y = F(z) can 
be inverted into z = F~ 1 (y). Accordingly, Eq. ( 113]) can be reformulated as 

exp[-/3?7 p (F- 1 (y))] 



n(y) 



-exp 



/3U(F~ 1 (y)) — /3U p (F~ 1 (y))\ , (15) 



for all ?/, < ?/ < 1. Eq. (|T5|) actually gives a definition of the function fl(y). The key 
result of H9] is that if 



n(y)~n y a -\ y^O, (a > 0), (16) 

the distribution of the sum z p = J2i=i l hi converges, up to a suitable rescaling, to one 
of the generalized extreme value distributions with parameter a. These distributions, 
illustrated in Fig. [2J belong to three different classes, depending on the large z behaviour 
of P(z). If P(z) decays faster than any power law (typically, a power-law potential U(z)), 
the generalized Gumbel distribution G a (x) is obtained, namely 

G a (x) = C g exp (— ay — ae~ y ^j , y — 6(x + v), (17) 

where 9 and v are rescaling factors introduced to have zero mean and unit variance, 
and C g is a normalization factor (see Appendix A). In the limit a — > oo, the generalized 
Gumbel distribution converges to the Gaussian distribution. If P(z) has a power law 
tail P{z) ~ z~ 1_fl when z — * oo, with /i > (corresponding, in the present mapping, to 
a logarithmic potential U(z)), one finds the generalized Frechet distribution, 

F *A%) = exp (~b f x^) , x > 0. (18) 

Finally, if P(z) vanishes above a certain value z max (say, there is a hard wall at z max ), 
and behaves as a power-law P(z) ~ (<2 ma x — 2;) M ~ 1 , for z — > z max , with ji > 0, then the 
resulting distribution is of the generalized Weibull type, 

Wa^x) = ax""- 1 exp {-b w x») , x > 0. (19) 

The parameters 6/- and 6 m are scale parameters, that can be tuned to obtain any of the 
distributions with unit mean or with unit variance (for the Frechet distribution, this is 
only possible if /i is large enough so that the mean or the variance are finite). 

As a rather generic illustration, let us consider the case when the potentials U(z) 
and Up(z p ) are given by 

U(z)=U z a , U p (z p ) = U^z;, (20) 

with a, 7 > 0. If a = 7 = 1, the particles and the piston are in a constant external force 
field, like the (local) gravity field, in which case Uq = mg and Uq = Mg. If 7 = 2, the 
piston is for instance linked to the hard wall situated at z = with a spring of stiffness 
k = 2U' . 

In order to determine the limit distribution of the volume fluctuations in the large iV 
limit, the parameter a characterizing the small y behaviour of £l(y) should be evaluated. 
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Figure 2. Examples of generalized extreme value distributions f(x), with a = 1.5. 
Left: Gumbel distribution f(x) = G a (x) (full line), Frechet distribution f(x) = F a ^(x) 
with /i = 1 (dashed line) and Weibull distribution f(x) = W a ^(x) with fj, = 1 (dot- 
dashed). Right: same distributions on a semi-logarithmic scale. 



Considering Eq. (|15[) . one needs to compute first the function F l (y) in the small y limit, 
which is deduced from the large z limit of F(z). Given that 

P(z) = \e~^ a , X= aiPU ( f ] /a , (21) 

with T(t) = Jq 00 duu l ~ x e~ u the Euler Gamma function, one has in the large z limit 

F(z)= rXe'^dx^ e -^ a (z ^ oo). (22) 

Jz r ( - ) z"- 1 



Inverting the relation y = F(z) to get z = F 1 (y), one finds to leading order in the 
limit y — ► 

F " lfe) K jmr- ( ln ? - 1 1 4) ln '% - ln r OT ■ (23) 

We now wish to compute f2(y). Putting Eq. fl23l) into Eq. f|T5|) . one needs to distinguish 

between the cases a = 7 and « 7^ 7. 
If a = 7, one gets 

l 1 ^) ( ^ 0) ' (24) 

with a = Uq/Uq. Therefore, Q(y) behaves as a power law in the small y limit, up to 
logarithmic corrections which do not modify the asymptotic distribution (see Appendix 
B). It follows that volume fluctuations are described by a Gumbel distribution with 
parameter a [19], which compares the relative intensity of the compression force acting 
on the piston and of the forces directly acting on the particles. In the limit of a strong 
compression force a — > 00, volume fluctuations asymptotically become Gaussian, as the 
generalized Gumbel distribution converges to the normal distribution in this limit. In 
the opposit limit a < 1 where the external compression force is small with respect to 
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the forces acting on the particles, the generalized Gumbel distribution converges to an 
exponential distribution (see Appendix A). 

In the case where the confining potential for the gas and that for the piston are of 
different functional form, a 7^ 7, one finds for Q(y), dropping logarithmic corrections as 
well as constants of order unity, 



Hence £l(y) does not behave as a power law when y — ► 0, so that the results of [19] do 
not apply. We show in Appendix C that when 7 > a, the limit distribution is a Gaussian 
law, while for 7 < a the limit distribution is exponential. This result is consistent with 
the following intuitive argument. When 7 > a, y£l(y) decays faster than any power law, 
and one expects this situation to be similar to the large a limit, for which a Gaussian 
distribution is recovered. In contrast, when 7 < a, yQ{y) decays slower than any power 
law, which is expected to be similar to the limit a — > 0, in which case one obtains an 
exponential distribution, as shown in Appendix A. The physical interpretation is that 
for 7 > a the piston confines the gas more strongly that the bulk confining potential, 
so that the latter becomes irrelevant at large size, leading to a regular confined gas with 
standard thermodynamic properties. For 7 < a the reverse is true; the gas is confined 
by the bulk potential; the piston becomes irrelevant and its fluctuations become those 
of a piece of flotsam driven by the fluctuations of the confined gas below. 

Note that we focused here on power law potentials, which correspond to a quite 
natural class of potentials. However, one could also consider logarithmic potentials U(z) 
and Up(zp), which would lead to Frechet distributions given in Eq. ( fT8l) for the volume 
fluctuations. Alternatively, if the potential diverges for a finite value z max , (for instance 
by adding a rigid wall on top of the piston), the asymptotic distribution of fluctuations 
would be of the Weibull type, as described in Eq. ( fl9l) . 

4. Discussion: relation with thermodynamics 

In the present paper, we have shown that the volume fluctuations of an ideal gas 
of classical and independent particles confined by an algebraic potential acting on 
both particles and piston along one dimension and by hard walls in the perpendicular 
directions, are described by generalized extreme value statistics. In the simple case 
when the piston is identical to the particles, the appearance of standard extreme value 
distributions is easily understood from a direct mapping of volume fluctuations onto 
an extreme value problem of independent and identically distributed (i.i.d.) random 
variables. 

In a more general situation the piston can differ from the particles in two ways; 
either the confining potential is of the same form but of different amplitude, meaning 
that the restoring forces on the particle and piston are different. Or, the functional form 
is different. In the first situation non-Gaussian height fluctuations still occur in the limit 





(25) 
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of large N in the form of generalized extreme value distributions parameterized by a 
real variable, a = U'q/Uq. If the restoring force on the piston becomes much larger than 
that on the particles, the distribution crosses over to Gaussian. Yet for fixed values of 
the forces, that is for fixed a, the distributions are non-Gaussian for all system sizes and 
no crossover occurs as a function of N. For different functional forms, a crossover does 
occur as a function of system size, either to Gaussian fluctuations if the piston is strongly 
confined, or to one-body non-Gaussian statistics, with exponential height fluctuations 
if the piston is less strongly confined. Our results are, in principle valid for a system of 
arbitrary scale perpendicular to the z axis. However, non-Gaussian fluctuations should 
be observable for a of order unity only, which implies a piston of microscopic extent 
perpendicular to z, ensuring that any experimental realization would be in the form 
of a quasi-one dimensional sample (another possibility could be to exert two different 
forces on the piston, that could be fine-tuned to compensate almost exactly). 

Non-Gaussian volume fluctuations have strong consequences for the thermodynam- 
ics. In previous work, non-Gaussian order parameter fluctuations have been related to 
critical phenomena or to the fact that an ordered phase is unstable in low dimensions 
[TT1 [28] . Analogous physics occurs for the models considered here: non-Gaussian vol- 
ume fluctuations lead to singular thermodynamics. To illustrate this we consider first 
the simplest case where piston and particles are identical, of mass m and confined by a 
gravitational force, mg, so that a = 7 = 1. In this case the height difference variables 
hj defined in (jSJ), are independent and exponentially distributed: 

P(hj) = (N-j+ l)f3mg exp[-(N - j + l^mghj], (26) 

as seen from Eq. (jHJ). One then has a 1/ f- noise- like spectrum [3], (hj) = ksT /(N — j + 
l)mg, giving directly 

(zp)«zoln(iV+l), ;?o = — , (27) 

mg 

where Zq is the characteristic length scale for the particles set by the gravitational field. 
Defining volume V = Sz p and external pressure, P = mg/S, leads to an equation of state 
for the confined ideal gas P(V) = kBTln(N) in the large TV limit. This singular non- 
extensive behavior signifies the crossover between a system confined by a bulk potential 
(the gravitational field on the particles) and an external constraint (the pressure imposed 
by the piston). It is mathematically equivalent to the case of a thermally excited one 
dimensional interface with long range interactions [3]. 

Let us now consider the case where the piston mass M is different from the mass 
m of the particles. In this case, the variables hj are still exponentially distributed, but 
now with 



mg(N — j + a) U m 

similarly to the "truncated 1 / /-noise" considered in [TS|[BJ]. The average piston position 
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(z p ) is then given by 

JV+l N+l 



(29) 



For large N, one can approximate the sum by an integral, yielding 

fN \ 
(z p ) « z In - + 1 . 



(30) 



Introducing again the (external) pressure P = Mg/S, one finds the equation of state 



where now a is a function of pressure. Hence, again we find that non-Gaussian 
fluctuations are associated with non-extensive thermodynamics. The ratio N/a = 
Nm/M compares the total mass Nm of the particles to the mass M of the piston (the 
mass M could also be an effective mass accounting for the constant force f p , positive or 
negative, exerted by an operator: M = M — f p /g)- If the piston becomes macroscopic 
with total mass, M, exceeding that of the gas, then we move into the regime where 
N/a <C 1. In this regime, we recover the ideal gas equation of state, P(V) = NksT, as 
well as Gaussian volume fluctuations. In the opposite case where N/a ^> 1 (typically 
if a is finite and N is large), one has non-Gaussian fluctuations as described in Sect. HJ 
and the non-extensive equation of state P(V) ~ akBTha.(N/a). Assuming that a is 
large but finite, and scaling the number of particles, N, with all other parameters held 
fixed, one therefore begins in the extensive regime for small (but macroscopic) N and 
V. In this regime (z p ) is less than zq and the effect of the confining field on the particles 
is negligible. On increasing N, one crosses over into the non-extensive regime when 
this length scale in exceeded. However, no crossover is observed in the statistics of 
fluctuations since a is large, and fluctuations are practically Gaussian even in the non- 
extensive regime. Taking q = 7 ^ 1 requires more calculation but leads to essentially 
equivalent results. 

A consequence of these results is that in the non-extensive regime the volume 
fluctuations are abnormally small on the scale set by the mean volume or the number 
of particles. Through the fluctuation dissipation relation, this scale is given by 



which, in the limit N/a ^> 1 is independent of (V). Hence the isothermal compressibility, 
k = — (l/(V))d(V)/dP, a normally intensive measure of the fluctuations, varies as 
l/ln(iV) and scales to zero in the limit, iV — > 00. Physically this result occurs because 
the potential confining the particles within the bulk of the sample suppresses the 
collective fluctuations present in the standard thermodynamic regime. The logarithmic 
dependence is characteristic of a marginal situation between the two regimes and is 
analogue to the marginal stability of an ordered phase at the lower critical dimension, 
such as the 2D-XY model [11] . or one dimensional interface with long range interactions 




(31) 




(32) 
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The case where the piston is more confined than the particles is best illustrated 
by removing the confining potential for the particles in the above example while 
keeping that for the piston. One now trivially finds the ideal gas equation of state, 
P{V) = NksT and regular thermodynamic fluctuations from Eq. Q. 



Appendix A. Large and small a limits of the generalized Gumbel 
distribution 



In this appendix, we wish to show that the generalized Gumbel distribution G a (x) 
converges to the Gaussian distribution when a — » oo, and to the exponential distribution 
when a — > 0. The distribution G a (x) is defined in Eq. (1171) . with 9, v and C g given by 



9 2 = *'(a) 



i[h« -*(«)], c g = ^_ 



The function \I/(a) is the digamma function defined as 

d 



ma) 



da 



lnlYa) 



(A.1) 



(A.2) 



where T(a) is the Euler Gamma function. 

Let us start with the case a — > oo, and determine the large a behaviour of the 
constants 9 and v given in flA.lj) . Using Stirling's approximation for the Gamma 
function, 



'27T 

T(a) « \l— a a e~ a 



oo 



one finds 



m:a) 

9 2 



-^-lnr(a) 

da 

1 1 

a 2a 2 



In a 



2a 



so that 9 and v are given for large a by 



6 = — 

>a 



0{a 



-3/2> 



2^d~ 

0(a-^ 2 ). 



0(a 



-3/2^ 



(A.3) 

(A.4) 
(A.5) 

(A.6) 



(A.7) 



It follows that 

9(x + v) = ^= + ±- 
y a 2a 

Hence for fixed x, 9(x + v) goes to zero when a — > oo, so that the term exp(— 9(x + v)) 
can be expanded to second order. Inserting the different asymptotic expansion given 
above in the generalized Gumbel distribution leads to, up to order a -1 / 2 corrections in 
the exponential, 



G a (x) 



exp 



If X 1 x l 

-\/ax a 1 = 1 

v 2 I ^ 2a 2a 



(A.8) 
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-x 2 /2 



2n 



[a — > oo . 



In the opposite limit a — > 0, the Gamma function behaves as T(a) 
^/(a) ~ — 1/a and 



a 



1/ = a In a + 1. 



The normalization factor C 9 reads 



„a In a 



1 



(a->0). 



9 r(a) 

Using 9{x + u) = (x + a In a + l)/a, one finds for G a (x) 



GJx) 



cxp 



fx + a In a + 1) — ae 



~ e exp 



_ e -(x+l)/a 



-{x+a In a+l)/a 



(a -> 0). 



(A.9) 
1/a, so that 

(A.10) 

(A.H) 

(A.12) 
(A.13) 



It is easily seen that the second factor in the r.h.s. converges to a Heaviside function for 
a -> 0, 



exp 



-(x+l)/a 



0(z + l) (a->0), 



(A.14) 



so that the generalized Gumbel distribution converges to the exponential distribution 
with zero mean and unit variance 



G a (x) -^e- (x+1) 6(x + l) (a->0). 



(A.15) 



Appendix B. Case = 7: Effect of logarithmic corrections 

The aim of this appendix is to show that the logarithmic corrections appearing in 
Eq. fl24|) do not change the asymptotic distribution. To that purpose, we follow closely 
the procedure of [TJ5], and keep essentially the same notations. Using Eq. ffTUl) and the 
results of [19], the distribution Qn(z p ) is given by 

K N 



Qn{z v 



-P(z p )Q(F(z p ))(l-F(z p )r. 



(B.1) 



-P' " AT! 

We assume for Q(y) a power-law form with logarithmic corrections when y — > 0, namely 



a-l 



In 



1 



(l/-0). 



(B.2) 



lnF(zp), we 



Let us define the value such that .F(z^) = a/N. In the limit N — > 00, diverges 
as -F(Zp) — > when z p — > 00. Introducing the auxiliary function g(z v 
perform the following change of variables, to look at fluctuations around z* N : 
* v 

Zn + 



Expanding g(z p ) in the neighbourhood of z* N leads to 

g( z p) = 9( z n) + v + e N(v), 



(B.3) 



(B.4) 
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where 

lim e N {v) = 0. (B.5) 

./V— >oo 



The distribution &n(v) is obtained from Qn( 

®n(v) = —^—Q n (z p ) 
9 \ z n) 

p) -F(z p )Q(F(z p ))(l-F(z p )f 



N\ g'(z* NJ 

where we have used P(z p ) = g(z p )F(z p ). In the large A limit, we have, keeping v fixed, 

F(z p ) « ± e -«, (B.7) 

(l-F(z p )) 7V ^exp(-ae-"), (B.8) 

g'(z p )/g'(z* N ) -> 1. (B.9) 

Using Eq. ( IB. 21) . one obtains for large A, as -F(z p ) — > 0, 

F(2p)n(F(«p)) « n (ln^ + 1;)' (B.10) 



~fi (^-j e— (lnA) d , (N-oo). (B.ll) 

Altogether, one finds 

M«) « §f(lniV) 5 fio (^)Y-exp(-ae^) (B.12) 



Let us now compute A^v, which is given by pj5] 

''dyn^Cl-y)^. (B.13) 



K N Jo 

With the change of variable v — u/N, we get for large A 
AM 1 /f, „/u\A uN* 1 



du ft — 1 « — / du — e _u . (B.14) 



it 



ATjv A jo VA/ V A/ A Jo VA 

Using the small y expansion of Q(y), one has 



A! 1 r N , ^ /mV" 1 A N\ a 
du iZo I — I I m 



A- A Jo VA/ V u 



~^( lniV ) 5r («), (tf-oo). (B.15) 
Coming back to the distribution $>n(v), one finally obtains from Eq. ( IB. 121) 

$n(v) -> $oc(^) = ^exp (-at; - ae"") , (A -> oo), (B.16) 
r(a) V ' 

which is precisely the generalized Gumbel distribution. In order to recover the standard 
expression G a (x) given in Eq. ([IT]) , one simply needs to introduce the normalized variable 
£ through v = 0(x + u), with # and i/ defined in Eq. (1A.1I) . 
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Appendix C. Case 0^7: Convergence of the volume distribution toward 
Gaussian and exponential laws 

In this appendix, we compute the asymptotic volume distribution in the case where 
£l(y) is given by 

(V-O), (C.l) 



£l(y) ~ — exp 
V 



with ?| = 7/a^l,?]>0 (we refer to Sect. [3] for notations), and b > 0. The derivation 
follows essentially the same steps as in Appendix B, but also bears some similarities 
with that done in Appendix A. The distribution Qn{z p ) is given by 



Qn(z p ) = 9'^) f(z p ) n(F(zp)) (1 - F( Zp )) N 



-pj — y \^pj - 1 \^p) ^v- 1 k^pjj v- 1 - - 1 k^pjj (C-2) 
Using the form (IC.lj) of Q(y) and recalling that F(z p ) = exp[— g(z p )], one finds for large 



Qn( z i 



K 



N 



g'(z p ) ex V [-bg(z p r] (l - 



We now make the following change of variables: 

* 1 



N 



V 



N ' br,g'(z* N )g(z* N )v-i 
where z* N is defined by g(z* N ) = lixN. One then has for large N 

-p 



The distribution $>m(v) reads 



K 



N 



exp 



-b(lixN) v -v - exp 



iV! 677(111 iV)^ 1 
In the case 77 < 1, (lnN) 1 " v — > 00 when A — > 00, so that 



briilnN)^ 1 



exp 



— exp 



brj 



e(v), (JV-oo), 



(C.3) 
(C.4) 

(C.5) 
(C.6) 
(C.7) 



and $tv(w) converges to the exponential distribution <&oo{v) = e~ v Q(v) (the prefactor 
converges to 1 by normalization of the distribution). 

In the opposite case rj > 1, the above argument does not apply anymore. Thus we 
start again from Eq. (1C.6|) . and make a saddle point calculation. Let us introduce the 
function iPn(v) such that $jv(v) = K N exp(— if) N (v)), namely 



4>n(v) = v + exp 



677 (In iV)''- 1 



with 



K 



K 



N 



N\br](\n N) 



^-.exphKlniV)^ 



We define v* through ij;' N (v*) = 0, leading to 



exp 



677 (In AO'' -1 



&77(lniV) 



n-x 



(C. 



(C.9) 



(CIO) 
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We then perform a change of variable 



v = v * + x\Jbr)(\n iV)^" 1 . 



(C.ll) 



For large N, we expand the term exp(— x/J 677 (In N) 11 x ) appearing in iPn(v) to second 
order, yielding 



then converges to a Gaussian law when N — > 00 (here again, the normalization of the 
distribution $at(o;) ensures that the prefactor converges to the correct limit). 
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The resulting distribution of x, 
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